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ABSTRACT 

We investigate models for the class of ultraluminous non-nuclear X-ray sources (ULXs) seen in a 
number of galaxies and probably associated with star-forming regions. Models where the X-ray emission 
is assumed to be isotropic run into several difficulties. In particular formation of sufficient numbers of 
the required ultramassive black-hole X-ray binaries is problematic, and the likely transient behaviour 
of the resulting systems is not in good accord with observation. The assumption of mild X-ray beaming 
suggests instead that ULXs may represent a shortlived but extremely common stage in the evolution of a 
wide class of X-ray binaries. The best candidate for this is the phase of thermal-timescale mass transfer 
inevitable in many intermediate and high-mass X-ray binaries. This in turn suggests a link with the 
Galactic microquasars. The short lifetimes of high mass X-ray binaries would explain the association 
of ULXs with episodes of star formation. These considerations still allow the possibility that individual 
ULXs may contain extremely massive black holes. 

Subject headings: accretion, accretion discs — binaries: close — X-rays: stars 
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1. INTRODUCTION 

The existence in spiral galaxies of off-nuclear X-ray 
sources whose luminosities appear significantly larger than 
the Eddington limit for a 1M Q object has been known 
for some time (Fabbiano, 1989). These sources are dis- 
tinct from the weak AGN present in many spiral galax- 
ies known as LINERs (Ho et al. 1997), although in at 
least one case (M33, cf Dubus et al, 1997) they may be 
confused with AGN. Recently considerable effort has been 
devoted to interpreting the properties of these 'ULXs' (= 
ultra-luminous compact X-ray sources, e.g Makishima et 
al, 2000 and references therein). A key to understanding 
their nature may be that they appear to occur preferen- 
tially, although not exclusively, in regions of star formation 
(Zezas, Georgantopoulos & Ward 1999; Roberts and War- 
wick 2000, Fabbiano, Zezas & Murray, 2001). In this paper 
we investigate models for the ULXs. 

Bright, non-nuclear X-ray point sources in galaxies divide 
into two groups: accreting neutron stars and black holes, 
and young supernova remnants. The luminosities of the 
first group, but not the second, are constrained by the 
Eddington limit: 



Edd 



AirGMinipC 



1.3 x 10 38 mi erg s" 1 , 



(1) 



where <tt is the Thomson cross-section and m\ is the ac- 
cretor mass M\ in M . This constraint applies to any non- 
explosive source, whether powered by accretion or other 
means such as nuclear burning. 

Evasions of the limit are possible, but rare. In the transient 
system A0538-66 a magnetic neutron star accretes from 
the wind of a Be-star companion. The system sometimes 
has super-Eddington luminosities Lx — 10 39 erg s _1 
(White & Carpenter, 1978) but these may result from the 
reduction in the electron scattering cross-section below 
ot in the ~ 10 11 G magnetic field of the neutron star per- 



vading the accretion columns. Thus if we regard the ob- 
served variability of ULXs as ruling out the identification 
as supernova remnants, a straightforward interpretation 
as non-explosive sources requires black holes with masses 
Mi ;> 50— IOOM0, accreting at rates which on occasion can 

exceed ~ 1O _6 M yr _1 . As we shall see, while individual 
ULXs could harbour such masses, there are major difficul- 
ties with such a picture as an explanation for the ULXs as 
a class. Accordingly we consider models in which the X- 
ray emission is assumed to be significantly beamed. In this 
case ULXs may correspond to a relatively shortlived but 
common epoch of the evolution of close intermediate- or 
high-mass X-ray binaries, perhaps the thermal-timescale 
mass transfer phase following the normal X-ray phase. 

2. LUMINOSITIES, LIFETIMES, MASSES, BIRTHRATES 

We first consider the restrictions placed by observation on 
accretion models for the ULXs. We assume that a com- 
pact object of mass M\ accretes from a mass reservoir 
(e.g. a companion star) of mass M%. We denote the mean 
observed number of ULXs per galaxy as n, the beaming 
factor as b (= fl/Air, where fi is the solid angle of emis- 
sion), the duty cycle (= time that the source is active as 
a fraction of its lifetime) as d, and define an 'acceptance 
rate' a as the ratio of mass accreted by Mi to that lost 
by M2, i.e. the mean accretion rate Mi = a (—Ma). We 
further define L S ph as the apparent X-ray (assumed bolo- 
metric) luminosity of a source, given by the assumption of 
isotropic emission, and let L40 = L sp h/10 40 erg s _1 . From 
these definitions it follows that the luminosity 



L = bL. 



P h 



10 40 &L 40 erg s" 



(2) 



and the minimum accretor mass if the source is not to 
exceed the Eddington limit is 



Mi > 10 2 fei 40 M G 



(3) 
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The total number of such sources per galaxy is 



(4) 



with a minimum mean accretion rate during active phases 
of 



Mi M 2 a , 

M activc = — i = J- > 10- 6 6L 40 M vr" 1 . 

a a 

The mass loss rate from M 2 is thus 

- Mi > 10" 6 — L ia M© yr" 1 , 
a 

and the lifetime of a source is 

M 2 6 m 2 a 

T = -M;^ 10 w^ yr ' 



(5) 



(6) 



(7) 



with m 2 
galaxy 



M 2 /M , leading to a required birthrate per 



40 



_ AT _n_ 6dL. 
~~ t" ^ M'10 6 m 2 a 



= 10" 



m 2 a 



40 _l 

yr . 



(8) 



It is important to note here that the required birthrate 
is independent of beaming (and duty cycle): the greater 
intrinsic source population N required by bd < ljci eq. 
is compensated by their longer lifetimes (cf eq. 



3. UNBEAMED MODELS 



For an unbeamed model we set 6=1, and recover from 
(pi) the requirement Mi > 10 2 L 4 o- We note immedi- 
ately that some ingenuity is required (cf Makishima et al., 
2000) to make these masses compatible with the charac- 
teristic observed X-ray temperatures 1-2 keV, whereas 
these are natural values for the lower masses we shall 
find below for beamed models. With 6 = 1 @ gives 
M ac tivc <; 10~ 6 £4o M yr" 1 . This effectively forces the 
mass reservoir to be a companion star; except for ex- 
tremely high black hole masses Mi £ 3 x 10 4 M Q , Bondi 
accretion from even a relatively dense part of the interstel- 
lar medium is inadequate, giving a rate 



M Bo „di a 1 x 10 



-11 "noo 



(p/m 



~ 24 g cm 3 ) 



Meyr- 1 (9) 
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where mioo = Mi/IOOM©, p is the mass density of the 
ambient interstellar medium, and uio,Cio are the rela- 
tive speed of the hole and the local ISM, and the local 
ISM sound speed, both in units of 10 km s _1 . While 
individual ULXs might contain black holes of masses 
Mi J> 3 x 10 4 Mq, it seems improbable that a galaxy like 
the Antennae should contain about 10 accreting exam- 
ples. Unbeamed models for the ULX class thus have to 
invoke a class of extremely massive X-ray binaries. As 
we shall see, this may be a potential problem. Accept- 
ing that suitable binaries could in principle form, there 
are the usual constraints familiar from LMXB evolution 
(cf Kalogera & Webbink, 1996). Particularly important 
are (a) the binary must be wide enough that the progen- 
itor of the compact star (here a ~ 100M Q black hole) 
must be able to fit inside its Roche lobe (otherwise it will 



provoke common-envelope evolution), and (b) the binary 
must be able to provide the inferred minimum accretion 
rate Mi ~ 1O _6 M yr -1 . Constraint (a) immediately sets 
a scale, as main-sequence stars of masses J> 100M Q have 
radii > 1O 3 R (e.g. Stothers & Chin, 1999). Using Ke- 
pler's law, and assuming Mi >> M 2 , this implies binary 
periods 

P Z l^Ioo 5 yr- (10) 

We can compare this with the critical period beyond which 
the accretion disc around the black hole cannot be ther- 
mally stable, and the system must therefore be transient. 
From King (2000) we find 



Pr, 



4TO 100 m 2 /S d - 



(11) 



We see that unbeamed ULXs must be transient. Hence 
the inferred Mi ~ 1O _6 M yr" 1 now refers to the out- 
burst state only: this is advantageous, as persistent mass 
transfer rates — M 2 of this order would have been diffi- 
cult to explain. To fill the Roche lobe in a binary with 
period given by (10) requires an extended star (note that 
the binary period essentially fixes the mean density of this 
star uniquely, cf e.g. Frank et al., 1992, Ch. 4). From 
the formulae of King (1988) (cf Ritter, 1999) we see that 
an evolved star with helium core mass M c ~ 0.4 — 0.5M Q 
will fill the Roche lobe, independently of the total donor 
mass M 2 > M c . Mass transfer is driven by the nuclear 
expansion of the star, at the rate 



M 2 ~ 1 x lO" 7 P y ° r 93 m^ 47 M yr 



(12) 



where P yr is the binary period in years (cf King, Kolb & 
Burderi, 1996, eq 7). Thus even a duty cycle d as long 
as 10% would yield mean accretion rates of the required 
order, i.e. Mi = M activc = Mid" 1 ~ 1O" 6 M yr" 1 during 
outbursts. 

However, applying the simple irradiated-disc theory of 
King & Ritter (1998) predicts that the accretion rate 
M ac tivc should decay linearly from its initial peak down to 
zero on a 10 -20 yr timescale. This is not easily compat- 
ible with comparison of ROSAT and Chandra data. This 
may not be a crucial objection to this type of unbeamed 
model for ULXs; the theory of outbursts in large irradiated 
discs is complicated, even without adding further difficul- 
ties such as radiation-induced disc warping. 
A more serious objection to unbeamed models is the one 
touched on above, namely that they require a ~ 1OOM 
black hole to coexist in a binary with an evolved compan- 
ion star. Mass loss from very massive stars with non-zero 
metallicity is usually thought to be so strong that the fi- 
nal black hole mass is much smaller than the initial stellar 
mass (see Baraffe, Heger & Woosley 2001 for a recent view, 
and Papaloizou 1973 for a possible objection). Assum- 
ing coeval formation of the two stars in the binary rules 
out a primordial origin for the black hole progenitor, and 
would thus require a progenitor with a mass ^> 1OOM . 
If the IMF within the stellar clusters is close to that of 
Salpeter, we conclude that the number of stars formed hav- 
ing masses > 1OOM is a factor ~ 100 lower than the num- 
ber of stars having initial masses 10 < m < 1OOM . We 
would therefore expect X-ray binaries containing neutron 
stars and lower-mass black holes to outnumber markedly 
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those containing higher-mass black holes. The X-ray lu- 
minosity of the systems observed in the Antennae (Fab- 
biano, Zezas & Murray, 2001) contradicts this, assuming 
that all systems are unbeamed, although the number of 
low luminosity systems is not currently well known. 
Alternatively, a ~ 100M Q black hole may only recently 
have gained a new stellar partner. Such black holes may 
be produced within dense clusters through the merger 
of lower-mass black holes (see for example Lee 1993, 
1995). Indeed this has been proposed as the origin of the 
moderate-mass black hole inferred to be present in the 
central regions of M82 (Matsushita et al 2000, Matsumoto 
et al. 2001, Kaaret et al., 2001). Although it is possible 
that a moderate-mass black hole produced in a central 
cluster of M82 has gained a stellar companion by some 
dynamical process (tidal capture or via an exchange en- 
counter involving a binary, this scenerio is unlikely to work 
for the systems observed in the Antennae where the ULXs 
are observed to be strongly associated with the young star 
clusters which are located some distance from the galac- 
tic nuclei (Fabbiano, Zezas & Murray, 2001). Any mas- 
sive BHs would therefore have to be formed within these 
stellar clusters and not the nuclear clusters. To produce 
moderate-mass black holes within a cluster via the succes- 
sive merger of lower-mass objects, the potential well of the 
cluster has to be sufficiently deep to retain the black holes. 
This can be the case for the stellar cluster in the nucleus 
of a galaxy but is not true for globular clusters where the 
typical escape speed is far too low to retain black hole bi- 
naries as they are hardened via encounters. This has been 
suggested as the reason for the absence of black hole bi- 
nary systems in globular clusters (Sigurdsson & Hernquist 
1993, Kulkarni et al 1993). 

A population of ~ 100M Q black holes originating from a 
much earlier generation of effectively zero-met allicity stars 
seems unlikely to explain the ULXs observed in the Anten- 
nae (Fabbiano, Zezas & Murray, 2001) as these black holes 
would be distributed throughout the galactic halo and the 
probability of picking up stars from the young stellar clus- 
ters via dynamical encounters within the last ~ 10 7 years 
is extremely low. 

4. BEAMED MODELS 

Since unbeamed models run into difficulties because of the 
required black hole mass ~ 100M Q and the need for a 
companion, we consider the effect of assuming that the 
observed X-rays are mildly beamed. The simplest candi- 
date mechanism is the idea that the accretion disc around 
an accreting black hole has much lower scattering opti- 
cal depth over a restricted range of solid angles (e.g. the 
rotational poles) than in other directions. Almost all the 
emitted X-rays would therefore emerge in these directions. 
A beaming factor b ^ 0.1 would bring the required mini- 
mum accretor mass (cf eqn ||) into the range M\ <J 10M Q 
commonly found in dynamical measurements of X-ray bi- 
naries, particularly quiescent soft X-ray transients (e.g. 
Charles, 1998), while b < 0.01 would bring Mi down to 
neutron-star values. In addition this kind of disc geome- 
try, i.e. a thick disc with a central funnel, may actually 
radiate a total luminosity in excess of the Eddington limit 
(Jaroszynski, Abramowicz & Paczynski, 1980; Abramow- 
icz, Calvani & Nobili, 1980). Thus such modest 6-values 



may allow quite large apparent luminosities for perfectly 
standard black-hole or neutron-star masses. The obvious 
implication is that beamed ULXs might represent some 
short-lived phase in the evolution of a large class of X-ray 
binaries: from eqn ^ we find r <^ 10 7 W2a/(6/0.1)c? yr. 
A good candidate for such a phase is an episode of 
thcrmal-timescale mass transfer. These are extremely 
common, occurring when the donor has a radiative en- 
velope and is either (a) somewhat more massive than the 
accretor, and/or (b) first fills its Roche lobe as it expands 
across the Hertzsprung gap. In general both cases give rise 
to highly super-Eddington mass transfer rates. Case (a) is 
unavoidable for example in any neutron-star binary with 
an intermediate-mass (~ 2 — 4M Q donor): King & Ritter 
(1999) and Podsiadlowski & Rappaport (2000) show that 
Cygnus X-2 is a survivor of such an episode, in which 
—Mi reached values of order ~ 10 _6 M Q yr" 1 and the ex- 
cess mass transfer is simply blown away from the system 
rather than resulting in common-envelope evolution (see 
also King & Begelman 1999, and Kolb et al., 2000). Case 
(b) requires only a reasonably wide binary separation after 
formation of the compact star, and clearly benefits from a 
large initial phase space. One of Case (a) or (b) is also the 
likely path for all high-mass X-ray binaries such as Cyg 
X-l once the current wind-fed X-ray phase ends. 
Until recently it has generally been assumed that thermal- 
timescale episodes are unobservable, as they are short, 
and without beaming X-rays could not emerge from the 
super-Eddington accretion flow at all. We investigate here 
the possibility that ULXs could be systems in this phase, 
where beaming allows us to see the X-rays. 
The thermal-timescale mass transfer rate from a donor 
near the upper main sequence is roughly (cf King & Begel- 
man, 1999) 

- m 2 ~ 3 x 10~ 8 m| 6 M yr" 1 . (13) 

Comparing with the Eddington accretion rate we can cal- 
culate an acceptance rate 

a = 0.43mim^ 2 - 6 (14) 

and thus a lifetime (assuming d = 1) 

t £ 4.3 x 10 6 (6/0.1)" 1 J L 40 1 mim 2 ' 1 - 6 yr (15) 

and birthrate 

B £ 2.3 x 10 _6 n(6/0.1)i 40 m]; 1 m^ 6 yr" 1 (16) 

per galaxy. In particular, for a system like Cyg X-2, which 
has mi ~ 1.4, m 2 ~ 3M Q (King & Ritter, 1999; Podsiad- 
lowski & Rappaport, 2000; Kolb et al., 2000) we get a 
required birthrate 

B^lx 10" 6 n(6/0.1)i 4 o yr -1 . (17) 

We may compare this with the Galactic birthrate ~ 10" 6 — 
10~ 7 yr -1 deduced for Cyg X-2-like systems (King & Rit- 
ter, 1999; Podsiadlowski & Rappaport, 2000; Kolb et al., 
2000). For a high-mass black-hole system like Cyg X-l 
both mi and m 2 are probably significantly higher, raising 
B by as much as an order of magnitude. However the short 
X-ray lifetime ~ 10 5 yr of this X-ray phase requires a cor- 
respondingly high Galactic birthrate ~ few x 10" 5 yr" 1 , 
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again allowing a significant ULX population. X-ray bi- 
naries reach the thermal-timescale phase in a timescale 
comparable with the main-sequence lifetime of the donor. 
Thus ULXs descending from high-mass X-ray binaries 
would naturally be associated with a young stellar pop- 
ulation, as required by observation. 

A possible example of a ULX in the Galaxy is 
GRS 1915+105, where L x ~ 1 X 10 39 erg s" 1 (e.g. Bel- 
loni et al., 1997). As this is a microquasar, with radio 
jet axis at about 70° to the line of sight (see Mirabel & 
Rodriguez, 1999), only mild beaming b ~ 0.6 is possible, 
even assuming that we view the system at the edge of the 
X-ray beam. However this is indeed sufficient to reduce 
the luminosity to sub-Eddington values. Moreover such 
a geometrical alignment is quite reasonable, as it offers 
an explanation for the very unusual long-term behaviour 
of GRS 1915+105. The system was not detected in X- 
rays until 1992, since when it has remained persistently 
bright with only short interruptions. The usual explana- 
tion of this as an accretion disc instability, prolonged by 
self-irradiation by X-rays (cf King & Ritter, 1998) would 
require an implausibly large disc mass. An attractive alter- 
native is that the X-ray light curve reflects slight changes 
in the X-ray beaming, which would have decreased enough 
in 1992 to allow us to see the X-rays. 

5. CONCLUSIONS 

We have considered models for the ULX class and reached 
the following conclusions. 

(i) Unbeamed models probably require a black hole of 
Mi J> IOOMq in ~ 1 yr binary orbit with an evolved donor 
star. Forming such a system presents considerable diffi- 
culties, and even then the likely transient behaviour of the 
accretion disc in such a wide system is hard to reconcile 
with observation. It is still possible that an individual ULX 
may contain a very massive black hole (Mi ;> 3 x 10 M ), 
perhaps accreting from the interstellar medium. 

(ii) The assumption of mild beaming (6 ~ 0.1 — 0.01) re- 
duces Mi to values already observed for Galactic X-ray 



binaries, and suggests that ULXs represent a shortlived 
phase of their evolution. The most likely candidate for this 
is the thermal-timescale mass transfer episode inevitable 
in a very wide class of intermediate- and high-mass X-ray 
binaries. This in turn suggests a link to the Galactic mi- 
croquasars (cf King, 1998, quoted in Mirabel & Rodriguez, 
1999). The short donor lifetime in high-mass X-ray bina- 
ries would explain why ULXs are associated with young 
stellar populations. 

The major theoretical uncertainty for (ii) above is whether 
beaming is a natural consequence of high accretion rates. 
Only large numerical simulations can address this ques- 
tion. Perhaps encouragingly for this type of model, not 
only are the X-ray spectra fairly similar to those of Galac- 
tic black hole systems, Kubota et al (2001) observed X-ray 
spectral transitions typical of such source in two ULXs. 
The same type of spectral and timing variability has also 
been seen in the X-9 source in the M81 field (La Parola et 
al [2001]). There are several possible observational tests 
of these ideas. First, continued X-ray monitoring with 
a view to detecting possible changes in beaming geome- 
try is clearly worthwhile. We note however that X-ray 
eclipses are unlikely in any beamed model, assuming that 
the X-ray beam axis is normal to the binary plane. Opti- 
cal identifications of ULXs might allow at least two kinds 
of test: if the total X-ray luminosities really are as large as 
predicted if there is no beaming, one might expect to de- 
tect photionization nebulae around ULXs. If on the other 
hand ULXs are beamed, and thus of normal stellar mass, 
one might hope ultimately to detect a spectroscopic period 
(say 10's of days) in a ULX within the Local Group. 
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